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ABSTRACT 

We examine the effect of self-gravity in a rotating thick-disk equiHbrium in the presence of a 
dipolar magnetic field. In the first part, we find a self-similar solution for non-self-gravitating 
disks. The solution that we have found shows that the pressure and density equilibrium pro- 
files are strongly modified by a self-consistent toroidal magnetic field. We introduce 3 dimen- 
sionless variables Cb, Cc, Ct that indicate the relative importance of toroidal component of 
magnetic field (Cb), centrifugal (Cc) and thermal (Ct) energy with respect to the gravitational 
potential energy of the central object. We study the effect of each of them on the structure of 
the disk. In the second part, we investigate the effect of self-gravity on the these disks; thus we 
introduce another dimensionless variable (Cg) that shows the importance of self-gravity. We 
find a self-similar solution for the equations of the system. Our solution shows that the struc- 
ture of the disk is modified by the self-gravitation of the disk, the magnetic field of the central 
object, and the azimuthal velocity of the gas disk. We find that self-gravity and magnetism 
from the central object can change the thickness and the shape of the disk. We show that as 
the effect of self-gravity increases the disk becomes thinner We also show that for different 
values of the star's magnetic field and of the disk's azimuthal velocity, the disk's shape and its 
density and pressure profiles are strongly modified. 
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1 INTRODUCTION 

The theory of accretion disks, motivated in a large degree by their 
occurrence in some binary systems, particularly cataclysmic vari- 
ables, has been most fully developed for the thin Keplerian disks 
(Pringle 1981 1. Based on their geometric shapes, accretion disks 
are generally divided into two distinct classes, thin disks and thick 
disks. The theory of thin accretion disks (Shakura & Sanyev 1973) 
is well understood whereas there is no universally accepted model 
for thick accretion disks. The current interest in this theory is due to 
the possibility that thick disks may be relevant to the understanding 
of central power sources in radio galaxies and quasars. Observa- 
tional evidence suggests that in the center of many galaxies, matter 
is somehow ejected to large distances and gives off high energy 
radiation as it interacts with the external medium. There are other 
theoretical reasons for pursuing the study of thick disks. In the the- 
ory of thin disks, radial pressure gradients are neglected and ver- 
tical pressure balance is solved separately.lt was shown that this 
approximation is valid as long as the disk is geometrically thin. 
This condition may be violated in the innermost region of accre- 
tion disks around stellar black holes and neutron stars. The study 
of thick disks provides a better theoretical understanding of thin 
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disks as a limiting case and enables us to deal with intermediate 
situations <FranketaLl992l 

When considering the formation process of astrophysical ob- 
jects, such as galaxies and stars, the most crucial factor is self- 
gravity. In the standard thin accretion disk model, the effect of 
self-gravity is neglected, and only pressure supports the vertical 
structure. By contrast, the theory of self-gravitating accretion disks, 
is less developed. Early numerical work of self gravitating ac- 
cretion disks began with N-body modelling (Cassen & Moosman 
1991 ; Tomley, Cassen & Steinman-Cameron 1991). The time evo- 
lution of non-self-gravitating viscose disk has long been studied, 
and now we have a good theory describing its steady structure 
and its basic time-dependent behavior ^Shakura & Sanyev 1973^ ; 
jLyndenbell & Pringle 1974} . But with the added assumption of 
self-gravity in the disk, it is not easy to follow its dynamical evolu- 
tion, mainly because the basic equations for the disks are highly 
nonlinear (e.g. Paczinsky 1978 1; I Fukue & Sakamoto 1992 1. To 
solve the nonlinear equations of self-gravitating disks, the tech- 
nique of self-similar analysis is sometimes useful. Several classes 
of self-similar solutions were known previously, but all of them 
considered a disk in a fixed, external potential. Self-similar behav- 
ior provides an important class of solutions to the self-gravitating 
fluid equations. On the one hand, many physical problems often 
attained self-similar limits for a wide range of initial conditions. 
On the other hand, the self-similar properties allows us to investi- 
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gate properties of solutions in arbitrary details, witliout any of the 
associated difficulties of numerical hydrodynamics. 

Pen (1994) presented a general classification of self-similar 
self-gravitating fluids. Fukue & Sakamoto (1992) also analyzed 
the vertical structure of self-gravitating disks, but it is impossi- 
ble to compare their models with realistic disks, because they 
computed the vertical structure using the thin-disk approximations 
for polytropes. Finally, using the numerically method proposed 
by Hachisu (1986) and Hachisu et al. (1987), Woodward et al. 
(1992) developed a code to investigate the interaction between a 
disk and its central object. However, they only considered poly- 
tropic disks. Hashimoto, Eriguchi & Muller (1995) presented a 
two-dimensional equilibrium model for self-gravitating Keplerian 
disks. They showed that the shape of the disk (or disk thickness) 
in flounced by the rotational law and the ratio of the disk mass 
to the mass of the central star. Bodo & Curir (1992)computed the 
equilibrium structure of a self-gravitating thick accretion disk by 
an iterative procedure which produced a final density distribution 
in equilibrium with the potential coming from it. They showed that 
the geometrical size and shape the disks influenced by self-gravity 
of the disk. 

Accretion disks, containing magnetic fields, have been 
the subject of intense study in recent years. The role of 
the magnetic field in the equilibrium of accretion disks has 
been investigated by some authors jslandford & zanj ek 1997^ ; 
jLub ow, Papalizou & Pringle 1994^ for the thin disk models , 
an ideal magneto hydrodynamics (MHD) equilibrium with az- 
imuthal velocity and poloidal magnetic field has been analyzed 
tLovelace et al. 19851 : jMobarry & Lovelace 1986) . Using numer- 
ical methods, these authors fotmd that the magnetic field may 
change the shape and angular momentum distribution in the 
disk. Thick disk configurations with a poloidal magnetic field 
has be en studied by ^Tripathy et al 19 9'0}, in the MHD frame- 
work fBanerj ee, Bhatt, Das, Prasanna; 1993} . They investigated 
the equilibrium structure of thick disks and their stability in the 
presence of a dipolar magnetic field due to a non-rotating central 
object. Their solution shows that the pressure and the density equi- 
librium profiles are strongly modified by a toroidal magnetic field, 
resulting from the interaction between the permanent dipolar mag- 
netic field and the inertia of the gas disk. In a magnetized disk, the 
inertia of the gas is expected to bend the magnetic field lines back- 
wards, creating a toroidal component, which in turn may collimate 
a hydrodynamic outflow over long distances, forming jets. They as- 
sumed that the disk is non-accreting, stationary, axisymmetric, non- 
viscous, magnetized and that it is in equilibrium around a compact 
object, with only an azimuthal motion V^. 

We are interested in analyzing the role of self-gravity in thick 
disk equilibrium in the presence of the dipolar magnetic field of 
a central star. The outline of this paper is as follows: the general 
formalism of the problem is discussed in section 2, a self-similar 
solution of the equilibrium of non-self-gravitating accretion disks 
in the presence of a dipolar magnetic field is constructed in section 
3, a self-similar solution of self-gravitating magnetized accretion 
disks is constructed in section 4 and a summary of the main ideas 
is given in sections. 



nore the influence of energy dissipation. We consider the disk as a 
non-accreting MHD flow around a magnetized , non rotating, cen- 
tral compact object. The disk is assumed to be stationary and axi- 
symmetric. We use a spherical polar, inertial, coordinate system, 
(r, 9, (j)), with origin fixed on the central object. The basic equa- 
tions are two components of the euler equation in (r, 8) direction 
and poison's equation: 
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where p, p, V^, B^ and denote the gas density, pressure, toroidal 
velocity component, toroidal magnetic field, and gravitational po- 
tential respectively. Also G and M are the gravitational constant and 
the mass of central star. Solution of equations (1-3), in general, is 
a difficult task. Therefore, to simplify the problem, we impose two 
constraints. We adopt for B^ and V0 the same form as that given 
in ( Banerjee , Bhatt & Das; 1995). They assume that the magnetic 
field of the central star is dipolar: 



Br = 2Bo(-fcos9, 
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where Bq is the magnetic field strength on the surface of the central 
star near the pole and R is its radius. They showed that the dipo- 
lar field with the azimuthal motion of the gas disk, could make a 
toroidal component of the magnetic field, that is: 



Ba, 
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where Bi is an arbitrary constant with the dimension of a magnetic 
field strength, and k ^ § is a real constant. The general solution 
for K/j is: 

^0 = ^o(^)<"+^'sin-2"0, (7) 
H 

where R denotes the radius of the star. Vo is a constant with the 
dimension of a velocity and n is a real constant. 

Now, we try to solve our equations with these constraints. At 
first we introduce x according to : 

r Rd 

where R and Rd are the star and the disk radiuses respectively. Now 
we rewrite the Eqs.(6)-(7) as: 
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where Bi,Vo are the strength of the toroidal component of the mag- 
netic field and the rotational velocity on the surface of the star re- 
spectively. We can see clearly that the structure of the magnetic 
field of the central star can be modified by a rotational velocity of 
the gas disk. 



2 GENERAL FORMALISM 

As stated in the introduction, we are interested in analyzing the 
role of self-gravity in a thick disk equilibrium in the presence of 
the dipolar magnetic field of a central star. For simplicity, we ig- 



2.1 MAGNETIC FIELD CONFIGURATION 

In this subsection we want to study the magnetic field of central 
stars in the presence of a rotating gas disk. In a magnetized disk. 
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Figure 1. 2D Magnetic field configurations in a X-Y plane as a function of 
9, which luns from 6min = to dmax = ti" ,at r = Tin = 10 and R=7 
for different value of /3 = 



Magnetic Field configuFation,^=0 Magnetic Field canflguration,fi=O.S 




Figure 2. 3D Magnetic field configurations as a function of 8, which runs 
from Bmin = to 9max = TT ,at r = r'i„ = 10 and R=7 for different 
value of = 



magnetic field lines are deformed by the rotating gas. In this case, 
the gas is expected to be tied to the magnetic field lines and its 
inertia causes them to be bent backward and create a toroidal com- 
ponent. 

To study the magnetic field configuration within the disk we 
will look at the magnetic field lines, which satisfy the following 
equation: 

dr rd6 

Br 
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In order to visualize the field line configurations, we now 
choose fc = — I in Eq.6. It is useful to express the results in a Carte- 
sian frame through the usual relations (X — r sin 6 cos 0, F = 
r sin 9 sin (j>,Z = r cos 9). If we apply these transformations to the 
force equation we can obtain the corresponding parametric equa- 
tions that generate the magnetic field configurations: 
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where "l>o is a constant of integration (because we consider axi- 
symmetry solutions, we can set it to zero without any loss of gen- 
erality), P ~ -g^ and rin is the disk inner radius. From the above 
equations we can see that the azimuthal component of the mag- 
netic field can affect magnetic field line configurations in the disk. 
In Fig(l) we plot several 2-D field line configurations. We can see 
that when the toroidal component of the magnetic field is larger, 
the deformation of the magnetic field line is more obvious. In Fig. 2 
we plot several 3-D magnetic field configurations with the same 
parameters. 

And now we come back to Eqs.(l-3) and introduce other di- 
mensionless variables for the gravitational potential, the density 
and the pressure of the disk's material: 
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By inserting these dimensionless parameters into Eqs.(l-2), 
and by performing some simplifications, we obtain: 
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where Cb, Cc, Ct are constants defined by: 
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By introducing three dimensionless parameters indicating the 
relative importance of the energy of the toroidal component of the 
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magnetic field (Cb), the centrifugal energy (Cc), the thermal en- 
ergy (Ct) with respect to the gravitational potential energy of the 
central object, we can study them separately. If this simplification 
is used for Poisson's equation, Eq.3, we find that: 



The Density profile 



The IDensity profile 
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where Md is the mass of the disk out to Rd- 

Mi = yi^dpo (25) 

We may introduce another dimensionless variable Cg that gives the 
importance of self-gravity. Cg is the ratio of the disk mass to the 
central object mass: 

3Md 
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We would like to investigate self-similar solutions that de- 
scribe the equilibrium structure of the thick disk. Thus we intro- 
duce self-similar solutions for the density and pressure of the gas 
disk: 
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The Iso— Density profile 



The Pressure profile 




(28) 



Since we do not consider energy transport mechanisms, no 
specific form for the equation of state is assumed. As a result of our 
model we find that the density and pressure clearly show no poly- 
tropic, isothermal or other simple relation. We show that the indices 
in the non-self-gravitating and self-gravitating cases are completely 
different and that the mass distribution of the disk changes due to 
this fact. Solving Eqs. (19-20) and Eqs. (27-28) gives e and 5. 



3 NON SELF-GRAVITATING SOLUTION 

In this section we try to derive a self-similar solution of non-self- 
gravitating disks. In this regime we put S equal to zero. If we put 
Eqs.(27-28) into Eqs.(I9-20) we obtain: 
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In order to obtain the indices of the self-similar solution, we 
require that in each equation the exponent of x be the same for all 
the terms. Here we get: 



2fc-4 + e = 2n + 2 = e 

And we have: 

e = 2 - 2fc 

n = -2 
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These equations state that for getting a physical solution, k 
must be less than unity. Then by putting the self-similar indices in 
the equations and by doing some simplifications, we obtain: 



Figure 3. To the top left: Variation of density along the meridional as a func- 



tion of 6 from 6„ 



: Oto( 



Trfor different value of Cc (rotational 



velocity). Top right: for deferent value of k, bottom left: Iso-density contour 
for different value of Cc, Bottom right: variation of pressure with different 
value of Cc 
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+ Cc sin** + (3 - 2fe)Ct — - 1 = 0(35) 
1 dP 



(2fc-l)Cs^^^^^^4^^+CcSin**ecote-Ct^^ = 0(36) 

We have two equations with two unknowns. The first equation 
can be used to express P as a function of D: 

^^^^ = ^33^^[(l-asin«e)I>(e) + (fc-i)CBsin-^'=0](37) 

With our definitions of po and Po, we have D{9 = |-) = 1 
and P{6 — ^) = \. If we introduce this boundary condition in 
Eq.37, Ct can be calculated: 



a=(^[l-Cc + (/c-i)C.], 

Now if we introduce Eq.37 into 36, we find: 
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The solution of this equation can be write as: 
Di(6i) = 3Cfl(2fc-l)/(e) + Doi, 
where: 



(45) 



I [I-Ccsm^ e']^ sm-^''-^ e cose de ■ (46) 
Jo 

Now we come back to our boundary conditions for the density. 
The maximum density is at ^ = where D{^) — 1. Now we can 
find the value of Dqi as: 

= [l-Ccsm^e]^\e=^D{^) = [1-Ce]^;(47) 

Thus: 

Doi = [1 - Cc]^ - 3Cs(2fc - 1)1(6 = |); (48) 
Hence we can find the density distribution: 

D{e) = [(1 - - 3Cfl(2fc - - 1(9))] 
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The important result from this equation is that the value of the 
density is not zero in the rotation axis. This value depends on our 
dimensionless parameters, Cb and Cc, where Cb represents the 
effect of the toroidal component of the magnetic field and Cc the 
effect of the rotation velocity. 

Fig(3) shows typical examples of the self-similar solution in 
equations (37) and (49). In Fig(3), the density profiles, the iso- 
density contours and the variations of pressure along the meridional 
axes are represented. The density profiles for different values of Cc 
(rotational velocity) and k are displayed. As can be seen, the pa- 
rameter Cc determines the overall shape of the matter distribution. 
By increasing the effect of the rotational velocity (Cc), the shape 
of the disk changes to a thinner one, which is due to centrifugal 
force. In Fig(3), we also present the iso-density contours for the 
same value of Cc- As can be seen, the shape of the disk is strongly 
modified for different values of Cc- We can see that in this case, the 
density does not vanish at 6 = 0, tt (in the rotational axes). Our so- 
lution agrees with jBanerjee, Bhatt, Das, Prasanna; 199?) solutions 
for such systems. 

We should mention that the solution we have found is a math- 
ematical one. As a result, the solution has a physical meaning only 
for parts of the parameter space. In this self-similar solution we 
have a singularity in 6 = 0,n that arises from sin 6 in the integral 
of equation (46); Thus we cannot study the effect of the magnetic 
field in the equilibrium structure of the gas disk. That is the na- 
ture of the solution. After we have found the analytical solution for 
non-self-gravitating disks, we will try to investigate the effect of 
self-gravity in such disks. 



4 SELF-GRAVITATING SOLUTION 

In order to gain a better understanding of the physical situation of 
the disk, we include self-gravity in the above equations. The self- 
similar technique is utilized once again. If we introduce Eq.(27-28) 
in to Eq. (19-20) we obtain: 
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For finding the self-similar solution we must set equal the 
terms that have the same power. We obtain: 
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Now we can identify the radial dependence of the density and 
pressure of the gas disk (e and 5 ) as functions of fc. We can insert 
Eqs. (50-51) into Eq.(24) with use Eq.(27) we can find: 
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Now if we introduce Eqs. (55-56) into Eq. (24) and do some 
simplifications, we obtain an ordinary differential equation with 
two unknown functions (D{9) and P{9))'- 
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In order to solve this equation, we need another relation be- 
tween D{e) and P{e)- By looking at the definitions of _ff & 
we can find: 
2 dH 



2fc - 1 de 



= w{e), 
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Now by solving Eqs. (59-60) we can find the angular depen- 
dence of the density and pressure of the gas disk (D, P). We lin- 
earize these equations, introduce another dimensionless constant 
defined as the ratio of the disk mass to the mass of the central ob- 
ject (Cg), and we use Q ~ instead of D in order to simplify the 
equations. Thus with some calculations, we get: 



dyi r,fo \ 

— = S[e,yi,y2,y3) 
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where J/i = Q = J/2 = P, J/3 = and S and T are complex 
functions of ^ , j/i, j/2, J/a. We solve these linear equations, with 
two point boundary conditions with the shooting method. Solving 
these equations gives the meridional component of the density and 
pressure. We used Naryan & Yi (1995) boundary conditions: 
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Cg=10, Cc=0.3, K=1/3 
CB=10 Solid line 
CB=1 Dashed lonre 
CB=0. 1 Dotted line 
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Iso— Density profile 
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Figure 4. Variation of the density of the disk for several values of the pa- 
rameters: the top one: is for different values of Cb that represents the im- 
portance of the magnetic field; the middle one: for different values of Cg 
that represents the importance of self-gravity; the bottom one: is for differ- 
ent values of Cc that represents the importance of the rotational velocity 



Figure 5. Iso-density contours of the self- gravitating disks: the top one: for 
different values of Cg that represents the importance of self-gravitation; the 
bottom: for different values of Cb that represents the importance of the 
magnetic energy. 
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By changing the value of the dimensionless parameters Cb, 
Cc, Ct and Cg, we can study the influence of the toroidal magnetic 
field, toroidal velocity, thermal energy and self-gravity of the gas 
disk on its own equilibrium structure. In Fig( 4-6) we plot some of 
our results. By comparing the results in Fig(4-5) which are calcu- 
lated for different values of Cg, we can see how the disk is more 
sensitive to the influence of self-gravity. Fig(4) shows the variations 
of the density for different values of Cg , Cb , Cc in meridional di- 
rection. Concerning the geometrical shape of the disk we can see 
the tendency of the disk to become thinner when self-gravity plays 
an important role (Cg is large). And then we can see that a strong 
magnetic field coming from the central star can produce a large 
toroidal component, which can change the thickness of the disk 
near the pole. Fig(4) shows the iso-density contours of the disk for 
different values of Cg and Cb ■ In Fig(6) we present the iso-density 
profile of the disk for different values of Cc and the iso-pressure 
profile of the disk for different values of Cb - We can see in the iso- 
pressure profile that near the pole there is a singularity; maybe this 
big pressure gradient could produce a jet from the pole! In order 
to test this effect we should investigate the time evolution of such 
systems. This could be the subject of future works. 

Concerning the ratio of the mass of the disk to the mass of 
the central star, Eriguchi & MuUer (1993) extended their code to 
include a central object for a thick disk and calculated the equilib- 
rium structure of the disk which rotates according to the j-const 
law, where j is the specific angular momentum. Comparing this 
with the toroidal star of Hashimito et al. (1993), where the central 
object is not included, the relative thickness of the disk is smaller, 
which is caused by effects due to gravity of the central star (see fig- 



ures in Eriguchi & Muller 1993). As for the rotation law, Hashimito 
et al. (1993) and Eriguchi & Muller (1993) obtained rather thick 
disks, while Hashimito et al. (1995) found a flat disk shape where 
the gravitational force of the central object becomes relatively week 
compared with the self-gravity of the disk. Our results agree with 
the latter one. We can see in Fig(4-6) that the shape of the disk (or 
disk thickness) is influenced by the rotational law and the ratio of 
the disk mass to the mass of the central star. 



5 DISCUSSION AND CONCLUSIONS 

In this paper, we find a global equilibrium of a self-gravitating disk 
around a magnetized compact object. We ignore the effects of en- 
ergy dissipation and radial flow. We have investigated a stationary 
model for a self-gravitating disk influenced by a strong magnetic 
field coming from the central object. We ignore the effects of the 
local magnetic field of the disk. We consider that the radial flow 
is negligible with respect to the azimuthal motion of the gas in the 
disk. We find two types of self-similar solutions for such disks. The 
first one is the non-self-gravitating self-similar solution that shows 
that the shape of the disk changes when the physical parameters are 
modified. The second one is the self-gravitating self-similar solu- 
tion. In this class of solutions we can study the effect of self-gravity 
on the shape and thickness of the disk. 

In the non-self-gravitating solution we find an analytical solu- 
tion for some parts of the parameter space. This solution is found 
with a self-similarity method, and it gives a picture of the equilib- 
rium structure of the gas material around a magnetized compact 
object. The solution shows that the rotational velocity of the gas 
disk can change the equilibrium picture. The important result that 
can be inferred from this solution is the non-vanishing value of the 
density on the rotation axis. This value depends on the dimension- 
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Figure 6. Iso-density contours of the self-gravitating accretion disks: the 
top one: for different values of Cc, which represents the importance of ro- 
tational velocity; the bottom one: iso-pressure contours of the disk for dif- 
ferent values of . which represents the importance of the magnetic field. 



less parameters, Cb and Cc, where Cb represents the effect of the 
magnetic field and Cc the effect of the rotational velocity. 

In the self-gravitating solution, we find the vertical (merid- 
ional) dependence of the gas density in the disk Fig(4). The exter- 
nal dipole field can change the thickness of the disk near the pole, 
and the existence of a sharp density gradient in the ^-direction indi- 
cates a rapid decrease of the matter density away from the equato- 
rial plane. But when self-gravity plays an important role in the disk 
equilibrium (by increasing Cg) we can see that the thickness of the 
disk decreases. 

In Fig(5-6) we show some iso-density and iso-pressure con- 
tours of the equilibrium configuration of the disk for different val- 
ues of Cb, Cc, Cg, k. In the iso-pressure profile we can see a large 
pressure gradient near the poles, which may lead to outflows in 
those regions! This idea would require further study, such as a dy- 
namical investigation. 

We see that the presence of self-gravity in a thick disk, can 
change the geometrical shape of the disk and plays an impor- 
tant role in equilibrium structure of the disk. We also see that the 
strength of the magnetic field can change the structure of the disk 
near the poles. Finally, we must admit that our model does not de- 
serve the name "accretion" disk, since we did not include the ac- 
cretion (mass) flow. 

This study, one of the first of its kind, is a search for an equilib- 
rium structure of a thick disk in the presence of an external stellar 
dipole along with a self-consistently generated toroidal magnetic 
field B^. The meridional structure of the disk is mainly due to the 
balance of plasma pressure gradient, magnetic force due to 5^, 
and centrifugal force. The existence of an equilibrium structure, in 
fact, encourages one to look for generalizations of the analysis to 
the case when the radial velocity Vr ^ 0, representing accretions. 
On the other hand, yet another important aspect to be considered 
is the generalization of the newtonian analysis to general relativis- 
tic formalism wherein space-time curvature produced by the strong 



gravitational field of the central object would modify the magnetic 
fields and introduce new features. 

Our disk model might also apply to systems where the central 
object is a black hole, a neutron star or an active galactic nucleus. 
However, in this paper we modelled the disk under the assumption 
of a Newtonian potential. For a system where the central object 
is a black hole, the pseudo-Newtonian potential will play an es- 
sential role near the central object. The modification arising from 
a pseudo-Newtonian potential to the disk will be investigated in a 
future paper Consequently, we cannot calculate the accretion lumi- 
nosity. Thus, to compare our results with other models or observa- 
tions, it will be necessary to include mass flow in our model. 

We are grateful to the referee for a very careful reading of the 
manuscript and for suggestions which help us improve the presen- 
tation of our results. We are grateful Mohsen Shadmehri for contin- 
uous encouragement and useful discussions. We thank Kattia Fer- 
riere, Daniel Reese and saeed farivar for their useful comments. 
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